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Abstract. The melting of a Wigner Crystal of electrons placed into a host polar material is examined as 
a function of the density and the temperature. When the coupling to the longitudinal optical modes of 
the host medium is turned on, the WC is progressively transformed into a polaronic Wigner Crystal. We 
estimate the critical density for crystal melting at zero temperature using the Lindeman criterion. We 
show that above a certain critical value of the Frohlich electron-phonon coupling, the melting towards a 
quantum liquid of polarons is not possible, and the insulator-to-metal transition is driven by the ionization 
of the polarons (polaron dissociation). The phase diagram at finite temperature is obtained by making 
use of the same Lindeman criterion. Results are also provided in the case of an anisotropic electron band 
mass, showing that the scenario of polaron dissociation can be relevant in anisotropic compounds such as 
the superconducting cuprates at rather moderate e-ph couplings. 
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1 Introduction 

As it is well-known, the formation of a single polaron is al- 
ready a many-body problem, since it deals with the inter- 
action of one electron with an infinite number of phonons. 
It is thus clear that treating a system of polarons at finite 
density is a very complicated task. Some theories have 
been proposed so far which neglect the Coulomb inter- 
actions, or treat them in the random phase approxima- 
tion. The very existence of polaronic bands ^ lies on this 
theoretical basis, although, to our knowledge, no rigor- 
ous justification of this approach has ever been given in 
the literature. Old but reliable results were obtained by 
Mahan for low electron-phonon (e-ph) coupling and high 
electron density y], and extended later to intermediate e- 
ph coupling by Lemmens et al. |^. Recently, De Filippis et 
al. M obtained interesting results valid at lower densities 
and higher e-ph coupling, giving a more complete physical 
insight to the problem in the metallic phase. 

A completely different approach of the same problem 
can be given starting from the low density regime, where 
the particles localize due to the long-range Coulomb re- 
pulsion [g||6|,|7|. In this limit, the treatment is somehow 
simplified because the statistical effects between the elec- 
trons can be neglected. The simplest way to introduce 
the problem we are dealing with is to consider a Wigner 
Crystal (WC) of electrons, which is the ground-state of the 
electron gas at low density, and dip it into a polar (ionic) 
host material. What are the effects of the polarization on 
the stability of the WC? 



A polar dielectric is a compound which is character- 
ized by two main sources of polarization: the electronic 
polarization, corresponding to the vibrations of the core 
electrons at optical frequencies (10^*~^^s~^), and the ionic 
polarization, which is carried by the longitudinal optical 



modes in the infra-red region {lolo 



10 
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), corre- 
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spending to the relative motion of the positive and nega- 
tive ions. An external charge moving in such material |q| 
interacts with the dipoles of the dielectric, and produces 
a polarization cloud which partially screens its Coulomb 
field. The composite state of an electron plus its induced 
polarization is called a polaron Q. Such a quasi-particle 
has a finite extension (the polaron radius), a finite self- 
energy (the energy required to form the bound state), 
and a renormalized mass Alp which is generally greater 
than the electron band mass. As customary and following 
Fi^ohlich |Q , we can define a dimensionless e-ph coupling 

as a = [m* /2h^ijjLo) e^/Sj which is the only relevant 
parameter for a single isolated polaron. Here m* is the 
electron band mass and the effective dielectric constant 
£ is defined as \/e — I/Eoo — l/ss, where Eg and Soo are 
respectively the static and high frequency dielectric con- 
stants of the medium. 

If we now consider two interacting polarons, two dif- 
ferent situations can occur. Under some very specific con- 
ditions, the two particles can attract and form a bound 
state, sharing the same polarization cloud. The proper- 
ties of such composite bosons — the bipolarons — have 
been extensively studied, especially in connection to their 
Bose condensation at low enough temperatures. However, 
bipolarons can only exist if the ratio -q = £00/ £s is lower 
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than a critical value rjc ~ 0.1, and for rather large e-ph 
couplings (a > 6 — 7) [0. On the other hand, for rj > 7]^, 
which is quite common in polar materials, two polarons 
repel at large distances as l/ssd. Throughout the present 
paper, we will be dealing with compounds where the net 
interaction between polarons is repulsive. 

As was stated previously, a reliable treatment of the 
many-polaron problem can be given starting at low den- 
sity, if we assume that the particles localize due to the 
long-range Coulomb interactions, forming a Polaronic Wig- 
ner Crystal (PWC). In that case, the polarization of the 
host material will react to screen the external charges, 
and this will change both their interaction and kinetic en- 
ergy, as compared to the ordinary WC of electrons. On 
one hand, the total Coulomb interaction energy is shifted 
from ^ e'^/Rs to ~ e^/ssRs, Rs being the mean distance 
between two electrons at a given density. On the other 
hand, the kinetic energy is also reduced owing to the po- 
laron mass, and behaves as ^ hj^ /AIpR^. From this simple 
argument, we see that the two effects above compete re- 
garding the WC stability. While the decrease of the inter- 
action energy favours delocalization, the decrease of the 
kinetic energy tends to stabilize the crystallized state. 

The correct evaluation of the balance between these 
two effects is thus essential in order to understand the 
melting mechanism occurring for increasing density. It is 
particularly instructive to consider the adiabatic limit ob- 
tained for vanishing phonon frequency lulq i<^ ^^ oo)- In 
that case, as given by standard polaron theory IQ , while 
the polaron radius and energy remain finite, the polaron 
mass becomes infinite, and there is no way to get mobile 
polarons. Therefore, at zero temperature, the PWC can 
only melt through ionization of the electrons from their 
polarization potential-wells. The latter must be screened 
out by the liberated carriers in the same fashion as for the 
usual Mott transition in non-polar semiconductors, and 
the system will eventually become metallic [Ol. We have 
already published a phenomenological Mott-like criterion 
which describes the insulator-to-metal transition (IMT) in 
this situation Hi: 



,1/3 (? 



(e/eoo) Rp ~ 0.25 



(1) 



where Rp is the polaron radius (the bound state radius), 
which is the relevant localization length in our problem. 

If one restores a finite phonon frequency {ujlo ¥" 0) as 
is the case in real materials, the problem becomes more 
intricate, since the polarons are a priori mobile (they have 
a finite mass). Upon increasing the density, the polarons 
themselves could be delocalized without being ionized, 
possibly leading to a polaron liquid beyond the transi- 
tion. The aim of the present work is to examine carefully 
the quantum melting of a PWC in all the e-ph coupling 
regimes, as was previously done for the more simple case 
of the electron crystal (which is recovered for a = 0). Some 
important results have already been published in several 
letters F 



151. We provide here the details concerning 
the model, approximations and calculations, and give a 
more complete description including results at finite tem- 
peratures and for anisotropic compounds |16|. 



The paper is organized as follows. In section 2, we 
derive a model for the insulating state which includes 
both clectron-phonon and electron-electron interactions, 
and point out the different approximations involved. Start- 
ing from the Frohlich model, we assume a crystallized state 
for the polarons, and expand the long-range Coulomb in- 
teractions between different particles up to second order 
in the small displacements relative to their equilibrium 
positions, neglecting anharmonicity and exchange effects. 
It can be demonstrated |lq,f6| that the resulting dipole- 
dipole interactions between different electrons are respon- 
sible for the dispersion of the collective modes of the PWC. 
Furthermore, above a certain critical density, the dipolar 
terms can lead to a phonon instability of the polaron lat- 
tice, corresponding to the spontaneous excitation of long- 
wavelength transverse modes. This phenomenon, as well 
as its consequences on the properties of the dielectric con- 
stant of the system, will be extensively analysed in a forth- 
coming paper ||I^ . 

In the present work, we will neglect such non-local 
interaction terms, restricting ourselves to the mean-field 
Wigner approximation. This corresponds to including the 
many-body effects through an effective local potential act- 
ing on each particle. Hence, the many-polaron problem is 
reduced to the problem of a single polaron in an external 
potential. Taking advantage of this approximation, in sec- 
tion 3 we study the PWC in the framework of Feynman's 
path-integral method [Q , which is known to give reliable 
analytical results at any values of the e-ph coupling. 

In section 4, we examine the insulator-to-metal tran- 
sition occurring for increasing density (at zero tempera- 
ture) in terms of the Lindemann criterion, as was done by 
Nozieres and Pines for the electron crystal [n9|. Since the 
polaron is a composite particle, we argue that there are 
two different Lindemann criteria to be considered: one for 
the center-of-mass motion, which describes the melting as 
in the ordinary WC, and one for the relative displacement 
between the electron and its polarization cloud, which de- 
scribes the ionization of the polarons. This allows us to 
distinguish between two differernt melting mechanisms, 
depending on the e-ph coupling strength. At weak or mod- 
erate couplings (a < a*), the situation is very much the 
same as in the bare electron crystal, i.e. the melting is 
driven by the increasing fluctuations of the localized parti- 
cles. On the other hand, at strong e-ph coupling {a > a*), 
i.e. at low but finite lulo, the crystal melting is due to the 
dissociation of the polarons. In that case, the electrons 
cannot carry with them their phonon clouds in the metal- 
lic phase, and the melting towards a degenerate polaron 
liquid can be excluded. We next extend the results to flnite 
temperatures, and derive an approximate T vs. n phase 
diagram for the many-polaron problem. For an applica- 
tion to the superconducting cuprates, we flnally generalize 
our treatment to systems where the electron band mass is 
an anisotropic quantity. The results presented in Table 1 
suggest that polaron dissociation could be a real physical 
possibility in such compounds. 
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2 Model and approximations 

In this section, we derive our basic model for the polaron 
crystal. Let us consider a polar material, characterized by 
the three parameters aforementioned: lolO: ^s and Soo- 
We next consider an ideal doping procedure, which con- 
sists in introducing into this "host" material a finite den- 
sity of electrons, which we express as n = (47ri?g/3)~^, 
plus a rigid jellium which exactly compensates the nega- 
tive charges. Owing to its polarizability, the host crystal 
responds to these excess charges. If we neglect magnetic 
effects, the total interaction energy is given by standard 
laws of electrostatics as lIT^ : 






E 



b d^r dt 



(2) 



where time integrations must be carried out along a path 
satisfying the equations of motion. If we suppose that the 
localization length scales of our problem are larger than 
the anion-cation distance in the host material, the latter 
can be treated as a continuous medium, and we can con- 
sider the total electric field E(r) and the electric displace- 
ment D(r) due to the doping charges (electrons plus jel- 
lium) as macroscopic fields. Such a treatment is standard 
in polaron theory, and leads to the well-known Landau- 
Frohlich model |§,0. The two quantities above defined 
are related to the polarization field P through the rela- 
tion D = E -I- 47rP, which can be used to eliminate E, 
and derive a Lagrangian which only depends on the ex- 
cess charges. Following Frohlich, we separate two sources 
of polarization: P = Pq -|- Pir- Po corresponds to high 
frequency oscillating charges and can be included through 
the optical dielectric constant £00- Pjr is due to the ionic 
distortion, and vibrates at a frequency ulo in the infra-red 
range. One obtains the following interaction Lagrangian 
(see ref. |1^ for more details about this derivation): 



L, 



1 



SttEc 



D^ 



D 



27r£ 



UJ 



LO 



P , >2 p2 



trons, one gets the following many-body Lagrangian: 



(3) 

where Pir and {A'ke/uj'J^q) Pir are conjugate moments, 
and integrations run over the entire space {d^r in the in- 
tegrands has been omitted). The different terms in eq. (ph 
represent respectively the interactions between the excess 
charges, between the charges and the longitudinal optical 
phonons, and the free-phonon Lagrangian. We can now 
specify the nature of the doping charges — say electrons 
and positive jellium, but the same applies to holes and 
negative jellium — by writing D= D^ + D, where 



D =EDr = -5: 



r-r,: 



(4) 



Yi being the electron coordinate. We will come back later 
on the precise formulation of the term D+ due to the rigid 
jellium. 

Separating explicitely the two parts of the electric dis- 
placement, and adding the kinetic energy of the free elec- 



^-E- 



27re 



'LO 



p2 2 p2 



(5) 



— /d+.d+ -^ /Vd-dt 

l-|D+.(ED-)+/p..-(EDr+D+) 



Aire 



The interacting terms in (pi) now correspond to the self- 
energy of the jellium, the Coulomb interactions among dif- 
ferent electrons, the electron-jellium, electron-phonon and 
jellium-phonon interactions. Since the jellium is assumed 
to be rigid, it is convenient to introduce the polarization 
field which only responds to the electron motion by mak- 
ing use of the following transformation: 



^ir — ^ ir . --'-' 



SO that (||) becomes: 



^-E' 



27re 



'LO 



Pir ~ ^LO^ir 



(6) 



(7) 



Jr/^^-^^-^^/Eor-D, 



SttEc 



Jv^J 



J^|D+.(EDr)+/p..-(EDn 



Aire 



(the phonon-jellium interaction has disappeared, and the 
remaining terms involving the jellium are now divided by 
£s)- The many-body Lagrangian in its final form (Q) de- 
pends on the following parameters: (i) the electron density 
n\ (ii) the parameters of the host material to*, Es, Eoo, ^lOi 
which determine the e-ph coupling a. 



2.1 Phonon integration 

The exact partition function of the system can be ex- 
pressed in the path integral formulation as: 



Z = Trie-I^") = / e^ V{path) 



(8) 



where path integrals run over all the electron coordinates 
and on the polarization field. 5* is defined as usual as the 
integral of L in imaginary time {t = it, (3~^ being the 
temperature, and taking h = fcs — 1). The trace oper- 
ation corresponds to integrating over all possible closed 
trajectories. Since the terms depending on the polariza- 
tion are quadratic, the field P^^ can be exactly integrated 
out iQ. As a result, the total effective action for the elec- 
trons is the sum of three terms: 

5({rJ) - SjMru^m + Se{{Vi}) + S,-ph{{v,}) (9) 
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The first term is constant at fixed density and does not 
depend on the state of the system: 



so that 



S. 



jellium 



(3 



Stte, 



D+ D^ 



(10) 



5e-p„({r,}) = Y. 






rP 



|i'»W-rj(o-)l 



drda 



(12) 



We have introduced the phonon propagator G^j (/3, t — a) = 
(n+l)e-"l^-'^l+ne"l^-'^l, together with n = (e'^'^-l)-^ 
Formula ( |l2|) represents the retarded interactions between 
electrons, mediated by the lattice polarization: the diag- 
onal terms {i = j) correspond to the interaction of each 
electron with itself, i.e. the polaron effect, while the off- 
diagonal terms {i ^ j) give a retarded attraction between 
electrons i and j. 

The above expressions (^-(^ are valid at any tem- 
perature. We will now apply our hypothesis of crystalliza- 
tion at low densities, and make use of several approxima- 
tions. 



2.2 Model for the Crystallized State 

Let us first consider the jellium. We will assume from now 
on that it is constitued of spheres of radius Rg with a 
uniform positive charge density p"*" = e/ (47ri?f/3). Each 
sphere, which carries a total charge -l-e, is supposed to 
be centered on the sites {Ri} of a Bravais lattice. This 
approximation is intermediate between considering actual 
doping ion potentials, and a jellium uniformly spread out 
in the system. There is a small overlap between adjacent 
spheres, that we will neglect in the evaluation of the elec- 
trostatic energy. This is convenient to evaluate the jellium- 
jellium and the electron-jellium interaction terms. In par- 
ticular, one can write: 



D+(r)=^D+(r) 



(13) 



where 



D+ (r) 



47ri?f 
e r R 



47r |r - Ri 



(r-R,) ;if |r-R,| < Rs 
; otherwise 



SttEc 



D,+ D^ 



3e^ 



bSsRs 



The second term in (|9|) is the electronic part: 

It contains the electron kinetic energy, the instantaneous 
Coulomb repulsion between the electrons, and the inter- 
action between the electrons and the jellium. 

The electron-phonon coupling effects are included in 
Se-ph, which is given by: 



8ne,^J ' J 2^£3R, -R, 



(14) 



We next assume that the electrons are localized around 
the sites Ri of the same Bravais lattice, and we introduce 
the small displacements u^ = r^ — Rj , so that 



1 



AwEs 



T>+ D ■ 



(»#j) 



£s Ri ~ R?" 



u, 



where we have defined cu^ = e'^/m*R^ — 07^/3, ujp being 
the electron plasma frequency. Both Se and Se-ph can be 
decomposed into diagonal (i ~ j) and off-diagonal (i 7^ j) 
terms, i.e.: 



= Se 



e 

'2Z 



■E 



«#j 



1 



|Ri -Rj +Ui(T) -Ui(T)| 



7dT 



E 



9e' 



lOe.R, 



-nUr) 



2 e. 



'-< (r) 



dT 



and 



Se—ph — / ^ 
i 



E- 



4£ io io |u»M-u,(a)| 



drda 



4£ Jo Jo |Rj -Rj +Uj(r)-Ui(CT)| 



-drda 



Our next approximation consists in expanding the ac- 
tion up to second order in the {ui}, thus ignoring all the 
anharmonic and higher order contributions. After some 
elementary algebra, one gets: 



5({u,}) = ys, + iy5. 



with 



S^ - -/?• 



9e^ 



13 



lOe.i?. JO 

2 rP rP 



2 ^^-^ "*■' 



2 2 £s 



Ga;i,o (/3, T - Cr) 



4e ./o Jo |u»('r) -u,(cr) 



drda 



and 



s.,^ 



Y,fA^;<unr)u]{T)dT 



(16) 

dT 
(17) 

(18) 



u^Loe\ J2 f f A-^G^,, (/3, T - a)uf (r) u] (a) drda 

a7 -^O -^O 



2e 
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where the indices a, 7 = {x, y, z) denote the cartesian co- 
ordinates. We have also defined the dipolar matrix ele- 
ments: 






i?5 



Expressions (|l6|)-(|l8|) constitute our basic model for the 
polaron crystal. Since we focus on the insulating state at 
low density, the action S can be treated semi-classically, 
i.e. neglecting the exchange between the fermions and 
quantum statistical effects. 



2.3 Comparison with the electron crystal 

Our model ( |lq ) was obtained by making use of the follow- 
ing approximations: (i) we consider a particular structure 
for the jellium, (ii) we expand the off-diagonal part of the 
action up to second order in the displacements u^ , and (iii) 
we neglect the exchange between different electrons. In or- 
der to analyse in more detail the consequences of these as- 
sumptions, we can compare our results with the previous 
treatments of the electron crystal. Carr pOl demonstrated 
that the energy per electron in a WC with bcc symetry can 

be expressed as an expansion in powers of Rs ■ Intro- 
ducing the dimensionless density parameter Tg = Rs/clq 
(oo = fi-^/TOeC^ « O.53A is the Bohr radius) and taking 
the Rydberg as the energy unit, the result reads 



E 



-1.792 2.66 



3/2 



- + 0{rf') + 



(20) 



with b < 1. Basically, our three approximations above con- 
sist in neglecting all the terms in powers of r^ lower than 
—3/2 in (|20|). In fact, the first assumption slightly over- 
estimates the Madelung energy, whose numerical value is 
lowered from —1.792/rs to — l.S/rg. The second hypothe- 
sis consists in neglecting the anharmonic {b/r'^) and higher 
order terms. These terms can become important close to 
the melting point, and can change its numerical deter- 
mination, but do not affect qualitatively the physics of 
the problem. Our last approximation is to neglect all the 
terms proportional to the overlap between the wavefunc- 
tions of different localized electrons. The corresponding 

energy terms in (20), i.e. the exchange terms, fall off ex- 

1/2 
ponentially with rs .As pointed out by Carr, they are 

negligible up to r^ ~ 10. In the polaron crystal, where 

localization is more efficient due to the e-ph coupling, it 

is quite reasonable to neglect them in the density range 

of interest (these terms are responsible for the magnetic 

properties of the crystal, that we do not consider here). 



3 The mean-field Wigner approximation: 
solution by the Feynman method 



is sensitive since, using Gauss' theorem, it can be easily 
shown that the total mean electric field of a particular 
sphere (averaged in time), including the jellium, the elec- 
tron, and the polarization field, vanishes for jr^ — R^j > Rs 
[pi. In other words, two different Wigner spheres can only 
interact through the deviations of Fj from equilibrium, 
and that goes beyond the mean-field approximation. Each 
sphere is then on the average an independent entity, and 
the effect of the other particles is already included in the 
harmonic potential of eq. (Eq). At this stage, all the elec- 
trons localized in their jellium spheres are equivalent and 
uncoupled, so that the system is fully described by any of 
the single-particle actions ( [17| ) — i.e. we can take S = St. 
As was stated in the introduction, the off-diagonal ele- 
ments Sij will be treated in a further publication. 

The retarded interaction in (|l7|), which is responsible 
for polaron formation, cannot be integrated exactly. Fol- 
lowing Feynman |lq] , we introduce a quadratic trial action 
of the form 



5*0 = — 



ii^dT — 



^w 



'dr 



^ I I \u{t)-u{s)\^GUP,r-<j)dTda, (21) 



whose parameters K and if will be determined variation- 
ally (here u^ has been replaced by the coordinate u of the 
only electron present). It can be seen that ( pi|) corresponds 
to the following two-body Lagrangian 

L, = !!!lu- + Ex^-^iu-Xf-^^u^ (22) 
° 2 2 2 ^ ^ 2 Es 

after the integration of X has been carried out. This model 
describes an electron subject to an external harmonic po- 
tential and bound with a spring to a fictitious particle of 
mass M and coordinate X. The latter represents the po- 
larization cloud associated to the electron (see Fig.|l]) . The 
constant K , which measures the strength of the retarded 
attraction, can be written a,s K — M'uP' . 



I 



/ 



electron 

* 
m , u 



/ / \ 



Wigner potential 



spnng 
K 



aux. particle 
M,X 



Fig. 1. The classical model corresponding to the action ([2l|): 
the polarization cloud is replaced by a fictitious particle of mass 
M. 



In the following, we will restrain ourselves to the diagonal 
part of the action S = ^^Si, neglecting the dipolar terms 
Sij. This is generally called the Wigner approximation. It 



The Lagrangian (g^) can be diagonalized into two har- 
monic modes Ui and U2, whose eigenfrequencies ai and 
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a 2 are given by: 



a 



1,2 



T 



(23) 
We have introduced the shortcut notation v ~ {K/ n) ' 
for the characteristic frequency of the two-body model 
(^ = rn*M/ (m* + M) is the reduced mass). The two 
eigenfrequencies obey the following relations: 

aia2 = wuw life's (24) 

Finally, the mass M of the auxiliary particle can be writ- 
ten as: 

,. / „,2X /^2 



m^ 



= 1-^ 



-1 



and the polaron mass is defined as: 

Mp = m* 4- M. 
If we introduce the quantities 

_ w' -oi{ _ ai-w" 

-ril n 9 5 ^2 o 2 ' 

the relation between the real coordinates and the normal 
modes is simply: 



U = ^lUi + ^2U2 



X = 



(ui -U2). 



(27) 



and the diagonalized Lagrangian can be expressed in canon- 
ical form: 

Lo - -Aiiil + -A2UI - -Aiajul - -A2alul (28) 



These two independent modes correspond respectively to 
the vibration of the polaron center-of-mass, and to the 
vibration of the electron inside the polarization potential 
well. We shall call them respectively the external and the 
internal degree of freedom. It can be seen from expressions 
( p9| ) that in the dilute regime, the external and internal 
energy scales are well separated, and the polaron can fol- 
low the vibrations imposed by the external potential as 
if it was a rigid particle. Basically, ai and a2 give the 
excitation spectrum of the polaron crystal in the frame- 
work of the Wigner approximation, as is sketched in FigJ3. 
(note that the electron can also be excited towards a free 



(25) 


overall polaron 

I(co) "^"'i"" 

1 1 


polaron internal 
transitions 


(26) 
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phonons 
(host lattice) 







Fig. 2. Schematic absorption spectrum of the polaron crystal 
within Wigner's approximation. The two energy scales ai and 
0:2 are well separated, and the (host) phonon frequency lies in 
between. 



state outside the polaron potential; the corresponding fre- 
quency, which is higher than a2, is not shown here). 



3.2 The crossover regime 



3.1 The low density regime 

At very low densities (r^ — > 00, luw ~ 0), the polarons 
are so far apart that their properties are unchanged from 
the single polaron case. In other words, the localizing po- 
tential V{u) = m*ijj^^u'^/2es acting on the electron is a 
perturbation with respect to the polaron energy, and the 
eigenfrequencies ( p3| ) can be expanded for small w^/Es, 
which gives: 



ai 



771* UJlY 



OJext 



02 -> V. (29) 

{iOext is defined as the frequency of vibration of a particle 
of mass Mp in a harmonic potential whose spring con- 
stant is m*uj'^/es, while v is the internal frequency of an 
isolated polaron pl[). The corresponding eigenmodes are: 



Ul 



U2 



R 



77l*U + MX 

'' Mp 
u X. 



(30) 



As was noted in the previous section, the localizing Wigner 
potential forces the polaron center-of-mass to vibrate with 
a frequency ujext, which increases with increasing den- 
sity. However, the polaron is a composite particle, and 
its velocity is physically limited by the characteristic fre- 
quency of the phonon cloud: when LOext approaches lo^o, 
the phonons are too slow to "follow" the overall motion, 
and only the electronic degree of freedom oscillates. In that 
case, both ai and 02 can deviate significantly from their 
asymptotic values (p9|). If the density is further increased, 
the excess kinetic energy is transferred to the internal de- 
gree of freedom and the concept of a rigid polaron vibrat- 
ing in an external potential breaks down. The crossover 
density is roughly given by the condition: 



that we express using (p9|) as: 



(31) 
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where rj — Soo/ss- In the limit of strong e-ph coupUng 
(a > 6), the mass of an isolated polaron is well approxi- 



mated by Mp/m* 
be expressed as: 



nie 



0.02a^ 18 



200- 



so that eq. (|3l|) can 



V 



ii-vY 



1/3 



(32) 



and the crossover density becomes independent on the e- 
ph coupling. As an example, if we take m* = nif., £s = 30 
and Eoo = 5, the crossover region corresponds to r^ ~ 20. 
In the opposite limit (a — *■ 0), the phonon frequency tends 
to infinity, so that oj^xt never reaches the value ujlo- In 
that case, the internal structure of the polarons can safely 
be neglected, and the polaron crystal tends to an ordinary 
electron crystal. 
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3.3 Results of the variational procedure 

As in the case of the single polaron problem, Feynman's 
variational method gives, for any a, the best analytical 
upper bound for the free energy F . This is obtained by 
minimizing the following expression: 



F = Fo--{S- 
P 

where the term Fq is defined as 



5n 



-/3Fo ^ I v{path)e- 



(33) 



(34) 



and (• • •) stands for / V{path) {. . .)e'^" / J V{path) e^" (path 
integrations must be carried out on closed trajectories). 
The calculation of the functional (p3) to be minimized 
follows the lines of Feynman's work (details are presented 
in the Appendix) . The ground-state energy E at zero tem- 
perature can be obtained by taking the limit /3 — > oo in 
the preceding expressions. The result is: 



E = Co + -{ai + a2 ~ w) 



3iw^ - al){al - w^) 



a 



4{ai + a2)w'^ 
dt 



^Jo y^|i-(l_e-.*) + |j(l-e-"^*) 



(35) 



where all frequencies are expressed in units of wloj a-nd 
Co — —{9e'^a'^ /5esrs){me/m*). The last integral is non 
elementary and must be calculated numerically, a 1^2 £^re 
taken as the two independent variational parameters, and 
w is related to ai^2 through (gj). The expression (|3|), af- 
ter minimization, represents the energy per particle of the 
PWC in the Wigner approximation. The result is sketched 
in Fig.H for different values of a. It is interesting to note 
that, upon increasing a, the evolution from a WC of elec- 
trons towards a PWC is gradual. In fact, all the curves 
in FigJ^ have the same behaviour, and the only differ- 
ence comes from the fact that the minimum of E becomes 
more pronounced and shifts to higher densities. This is 
not surprising, since it is known that the formation of 
large polarons is not a real phase transition, but rather a 
continuous crossover 1231. 



Fig. 3. The energy per polaron in the crystallized state versus 
Ts (the single polaron value Ep, corresponding to the limit 
rs — 00, has been subtracted out). From top to bottom a = 
(electron crystal), a = 3, 5, 7, 10. The energy unit is meV. In 
this and the following figures, the parameters of the host lattice 
are Sa = 30, Eoo = 5, m* = me. 



3.3.1 Eigenfrequencies 

The three frequencies w, ai and a2 are illustrated in Fig.H 
as a function of the density, for two different values of a. 
As noted above, at weak e-ph coupling, the phonon cloud 
is fast and it can easily follow the polaron motion at any 
density. The external frequency ai is proportional to uw, 
while a2, which tends to w in this limit, is almost density 
independent (see inset). The two frequencies are always 
well separated, and the polaron can vibrate as if it was a 
rigid particle. 

On the other hand, at intermediate and strong q, the 
internal and external degrees of freedom are fully decou- 
pled only at low density (see main plot). When we reach 
the crossover region (indicated by an arrow), the exter- 
nal frequency "saturates" around wlo'- the phonon cloud 
cannot follow the vibrations imposed by the Wigner po- 
tential, and the polaron center-of-mass is virtually frozen. 
Upon increasing the density, the excess kinetic energy is 
thus transferred to the internal degree of freedom, and a2 
becomes proportional to low. 



3.3.2 Polaron mass 

In Fig.pl we illustrate the behaviour of Alp as a function 
of Ts, for different values of a. It can be seen that Mp 
always increases with the density when one approaches 
the crossover regime. A similar effect is well-known for 
an isolated polaron moving at a finite velocity v [p3| . 
If the quasi-particle kinetic energy Mpv^/2 reaches the 
phonon energy, the effective polaron mass increases be- 
cause phonons can be emitted from the polarization cloud 
in an incoherent way (the equivalent of Cherenkov effect in 
transparent media) . If we consider a polaron in a Wigner 
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Fig . 4. The characteristic frequencies of the two-body model 
(E2), in units of uilo- ct2 (upper curve), w (dashed curve), 
ai (bottom curve). The e-ph coupling is a — 10. The arrow 
marks the crossover region (see sec. ^j). Inset: same as main 
plot, with a — 1. 
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Fig. 5. The polaron mass Alp, normalized by the single po- 
laron value, for a — 3,6, 10. 



sphere, whose average velocity is v « ujexty^jMpuj^, it 
is easy to verify that the condition Mpv^/2 « huij^o gives 
exactly the crossover condition (pl|). 



4 The insulator-to-metal transition 

The melting of the ordinary electron crystal at zero tem- 
perature has been extensively studied in the literature, 
and there exist many different criteria which allow to eval- 
uate the critical density for the transition (see Care |24| for 
a review). A particularly simple approach was applied by 
Nozieres and Pines |1^ , which is based on the Lindcmann 
criterion. This phenomenological criterion, generally used 



to describe the thermal melting of the atomic lattices, 
states that the crystal is unstable when the spatial fluc- 
tuations of each particle around its equilibrium position 
exceed some fraction of the inter-particle distance. Such 
criterion will be generalized here to the case of a PWC, 
and will be used at both zero and finite temperature in 
order to describe qualitatively the insulator-to-metal tran- 
sition. 



4.1 Lindemann Criterion 

Calling u = {x, y, z) the coordinate of a localized electron 
in an ordinary Wigner crystal, the Lindemann criterion 
states that the crystal melts when 

{Suy/^/Rs > S 

where (5 is a phenomenological constant usually taken as 
S ~ 0.25. In the Wigner approximation, the average zero- 
point displacement in each space direction (say x) is given 
at zero temperature by: 



{6x^) 



2\l/2 



2meUJ\Y 



1/2 



(36) 



If we consider three dispersionless modes of equal fre- 
quency LOW = (^p/y/o, as is the case if we neglect dipole- 
dipole interactions, we get a critical density parameter 



.w _ 



64, which is close to the Monte-Carlo result (r, 



w _ 



100±20 in 3D, 123]). Of course, this criterion is only quali- 
tative and the resulting melting density is very sensitive to 
the choice of d, and moreover it does not give any precise 
information about the nature of the ground-state beyond 
the melting point. However, especially in the polaronic 
case, we shall see that it gives a good insight in the physics 
of the problem. 

The treatment that we have described in the previ- 
ous sections, which reduces the many-polaron problem to 
the problem of a single polaron in an effective potential, 
gives reliable results for the insulating state at any value 
of the e-ph coupling. It was shown in section that the 
system passes continuously from an ordinary Wigner lat- 
tice at weak a to a lattice of polarons at strong a, with no 
symmetry breaking between the two limits. However, we 
will show here that a definite distinction between the two 
cases comes from the melting mechanisms. Indeed, tak- 
ing into account the composite nature of the polarons, we 
can introduce two different Lindemann criteria, whether 
we analyse the fluctuation of the polaron as a whole, or 
rather the fluctuation of the electron with respect to the 
polaritazion field. Using the definitions ( po| ) of the center- 
of-mass R and of the relative coordinate r, we can write 
the following conditions: 

(i) {SRY^^/Rs> 0.25 
(u) {SrY^yRs> 0.25 

The former gives the critical value re for the melting 
towards a polaron liquid. The latter gives the critical value 
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Tc for polaron dissociation. The calculation of ((5i?^) and 
(fr^) in terms of the variational parameters is presented 
in the appendix. The result at T = is (all frequencies 
are expressed in units of uj^q)'- 



{SrY^'/Rs = 



When the e-ph coupling vanishes {Mp — > to* and R — > 
u), eq. (i) reduces to the ordinary Lindemann criterion 
for the WC, and the transition is due to fluctuations of 
the (weakly renormalized) electrons. On the other hand, 
the fluctuation of the internal coordinate gives a measure 
of the polaron radius [ p6[ , so that criterion (m) can be 
written as Rp/Rs = const, which closely resembles eq. 
(H) above. As a matter of fact, we can consider the second 
criterion as a generalization to finite phonon frequencies of 
the dissociation argument given in the introduction (which 
we have shown to be valid in the limit lolo ~^ 0, when the 
potential acting on the electron is completely static). 



4.2 Transition at T=0 

We are now in a position to analyse the IMT for any val- 
ues of a. Fig.ra shows the calculated ratios which enter in 



a) 1021 



n (cm-3) 
lO'" lO'' 



10'; 



n (cm-3) 
10i» 1017 1015 





.a=l 


0.8 


\ 


^0.6 


a=3\ 


A 

io 0.4 
V 


-?^^ 


0.2 


-a=10 -^^^ ^ 


n 


a=100 



10 



100 
r. 



1000 




1000 



Fig. 6. The ratios which enter in the Lindemann criteria 
(i) and (ii): a) {SRY^'^/Rs for a = 1,3,5,7,10,100; b) 



(5. 



,2\l/2 



/Rs for a — 1 (dashed line), a = 3 (dotted line), 



Q = 10 (full line). The horizontal line corresponds to the phe- 
nomenological melting value S = 0.25. 



the conditions (i) and (ii), as functions of r^, for differ- 
ent values of a. One observes that for a fixed a, while 
{Sr^y^'^ /Rs always increases with increasing density (all 
the curves in FigMb have roughly the same behaviour), 
the center-of-mass fiuctuation {SR'^)^^^ /Rs becomes bell 
shaped for sufficiently large a (it changes its slope exactly 



where, according to the argument of section 3.2, the in- 
ternal structure of the polaron becomes important — see 
Fig.^.a). In addition, the maximum value of {6R'^)^/'^ /Rg 
decreases with increasing a. Remembering that the tran- 
sition occurs when one of the fluctuation ratios reaches 
the value 0.25, we see that there exists a certain criti- 
cal value a* above which the crystal melting a la Wigner 
is prevented, and the driving mechanism switches to the 
polaron dissociation, typical of the static limit. An impor- 
tant consequence of this result is that, for a > a*, the 
polarons do not survive beyond the IMT. In other words, 
in the strong e-ph coupling limit, a liquid state of large 
polarons cannot exist at zero temperature. With the pa- 
rameters of Fig. pi we find a* ~ 7.5. 
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Fig. 7. Zero temperature phase diagram as deduced from cri- 
teria (i) and {ii): Tc for crystal melting (full line) and r" 
for polaron dissociation in 3D (dashed line). The correspond- 
ing quantities for an anisotropic polaron crystal (see section 
4.4 ) are drawn in grey. The dotted line corresponds to the low 
density approximation (BSl). The melting driven by the fluc- 
tuations of the localized particles (criterion (i)) is impossible 
beyond a* (full circle). The critical density given by eq. (hi) is 
indicated by an arrow. For regions I, II and III, see text. 

Fig.^ illustrates the approximate phase diagram ob- 
tained from the two Lindemann criteria at T = 0. The 
full line corresponds to criterion [i). We see that r)/ de- 
creases with increasing a (the corresponding critical den- 
sity increases). This can be easily unde rstood in terms of 
the asymptotic behaviour of section 3.1, As we have stated 



above, at low density, we can neglect the composite na- 
ture of the polaron and its mean square displacement in 
an external potential of frequency ujpxt is given by: 



{5Ri- 



2\l/2 



2MpL0e 



(37) 



which fits very well the low density part of {SR^Y'"^ /Rg 
(right side of Fig. |[(a)). In the weak coupling regime (in 
practice for all a < a*), we can use this formula to cal- 
culate the critical r^ for crystal melting. The result can 
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be expressed in terms of the critical parameter r^ of the 
ordinary electron crystal as: 



M 



w 






(38) 



The competition between the kinetic and potential energy 
is clearly visible in eq. (^8|) : if the mass of the carriers is 
strong enough, despite the static screening of the charges 
which reduces the Coulomb interactions, the crystalliza- 
tion is favored with respect to the ordinary WC of elec- 
trons. The melting curve corresponding to this approxi- 
mation is shown in Fig.0 (dotted line: we see that ( p8|) 
slightly underestimates the critical density, especially for 
a close to a*). Since the polarons do not dissociate at 

the transition, we expect the metallic state beyond ric 
to be either a renormalized electron liquid (region I) or a 
degenerate polaron liquid (region II). 

For a > a*, the IMT is described by the dissociation 
mechanism (ii) , and the critical density is almost constant 
(in a nutshell, the crystal is destroyed when the polarons 
start to overlap). In this range, although the polarons do 
not survive beyond the transition, the resulting electron 
liquid should present both strong e-ph and e-e interactions 
(region III) . A second possibility is that not all of the po- 
larons are destroyed at the transition, because this would 
cost too much potential energy, especially in the strong 
coupling limit. This can be understood as follows: when 
the density of localized polarons is increased, there is a 
competition between the increase of {Sr"^)^^"^ /Rs due to 
the zero-point fluctuations (which destabilizes the PWC), 
and the tendency of the system to preserve the poten- 
tial energy stored in the polaron bound states (which in 
turn should prevent dissociation). A compromise could be 
achieved in a hypothetical mixed phase, where some itin- 
erant carriers could coexist with the localized polarons. In 
fact, some experiments p^ suggest that all the carriers 
introduced beyond the critical density are itinerant (their 

number being proportional to n — n" ), the density of 
localized polarons being kept constant and equal to Uc ■ 
The properties of such a mixed phase are currently un- 
der studies. It is clear, however, that this could be stable 
only close to ric , while for n ^ Uc , all the interactions 
would be screened, and the system would tend to a normal 
Fermi liquid. 



4.3 Transition at T ^ 

If we go to finite temperatures, the melting of the polaron 
crystal can be studied in terms of the same Lindemann 
criteria (i) and (m) as in the zero temperature case, pro- 
vided that we include the effect of thermal fluctuations 
in the definition of ((5i?^) and {6r'^). The correct result 
is given in the appendix (see eq. IJI(\ ) and (|7l|)), where 
the variational parameters ai and a2 are now obtained 
by minimizing the free energy F given by eq. (Eq), to- 
gether with (pl|), ( |62| ) and (|6q). It is known from the sin- 
gle polaron case, however, that the upper bound (p3|) that 
one estimates starting from the two-body model (63) only 



gives correct results for T < lulq P^^M- ^t higher tem- 
peratures, such quantities as the polaron mass Mp and 
internal frequency a2 are ill-defined, and the treatment 
is no longer applicable in this simple form. Nevertheless, 
since the characteristic energy scale for the thermal fluc- 
tuations of the polaron center-of-mass is ai < lulq, we 
reasonably expect the thermal melting of the crystal to 
occur at temperatures well beyond the limits of validity 
of the treatment, as can be verified a posteriori. 

In order to understand the basic physics of the prob- 
lem, we first restrict ours elves to the low-density regime 
as was done in section 3.1. This allows us to express {SR'^) 
and ((5r^) by the approximate formulas 



((5i?2) = (2Mpai)-icoth^ 
{6r^) ^i2iia2r'coth^. 



(39) 
(40) 



Polaron theory predicts that at temperatures much lower 
than ujlo i the properties of each polaron are almost unaf- 



fected by thermal fiuctuations 1 28 



In that case, we can 
consider to a first approximation a2{T) w q;2, Mp(T) 
Mp, and consequently ai{T) « a-t 



Yi (fromeq. (g)). There- 
fore, in this simple limit, the only temperature dependence 
of the radii (^9|) and ( [40|) comes from the explicit /3 factor. 
In the intermediate and strong coupling regimes (i.e. 
where the mechanism of polaron dissociation can become 
relevant), since a^ ^ 'jJlo^ coth(/3a2/2) k. 1, and we im- 
mediately see that the dissociation density does not de- 
pend on the temperature: 



„('*) 



(T)«rf)(T = 0) 



(41) 



On the other hand, criterion (i) can be solved to give the 
critical temperature for crystal melting: 



Tc = Wp. 



At low temperatures, this equation has two different solu- 
tions which correspond respectively to the quantum melt- 
ing analysed in section 4.1, at r^ « re {T — 0), and to the 
classical melting due to thermal fiuctuations, at rg — > cx). 
Let us first expand eq. (^) in the classical case: 




2% .2 1/3 



e.^r. 



(43) 



The melting temperature in that case is proportional to 
the only relevant energy scale in the problem, i.e. the 
strength of the Coulomb interactions [p9[ . In the quan- 
tum limit, i.e. for r^ close to the zero temperature value 



„(») 



(T = 0), Tc can be expressed as: 



T, 



-Wext/ln 



Irs 
4 



(i) 
~ rl'{T=o) 

re (T=o) 



(44) 



which states that the critical temperature tends to zero at 



Si) 



(T = 0) with a slope that diverges logarithmically. 
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In the general case (i.e. at any density and any e-ph 
coupling), the Lindemann criteria (i) and (ii) can be ac- 
curately evaluated as was indicated at the beginning of 
this section. In figures || and ||, we show two characteristic 
phase diagrams in the (T, n) plane, which correspond re- 
spectively to e-ph coupling constants above and below a* . 

(Fig.|, 



In the strong coupling regime, where ric > 
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Fig. 8. The T vs. n phase diagram for a = 10 > a* 
obtained from criteria (i) and (ii). 
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Fig. 9. The T vs. n phase diagram for a = 7.2 < a*, as ob- 
tained from criteria (z) and (ii). The shaded area corresponds 
to temperatures T < Tdeg, where the polaron liquid could be- 
come degenerate. 



diagram is divided in three regions, but now the polaron 
liquid can become degenerate at sufficiently low tempera- 
ture (i.e. for T < Tdeg), in the intermediate region between 

(shaded area in Fig.O). 



(i) 1 (ii) 

ric and Uc 



a = 10), while quantum fluctuations lead to polaron disso- 
ciation, thermal fluctuations lead to the melting through 
criterion (i), so that a polaron liquid state can be achieved 
by increasing the temperature above Tc- However, a simple 
argument shows that it is very unlikely that such polaron 
liquid is degenerate. In fact, the particle statistics become 
relevant when the de Broglie thermal length is comparable 
to their average distance, i.e.: 



2^2 



MpksT 



1/2 



2i?,, 



(45) 



Therefore, the polaron liquid is expected to behave clas- 
sically down to temperatures of the order 



2 • lO^K 



Ideg 



Mf 



(46) 



4.4 Extension to the anisotropic case 

Both the static (Es) and high frequency (Eoo) dielectric 
constants have been measured in some undopped cuprates. 
For the parent compound La2Cu04, it was found [ pO| 
Eg ~ 30 and £oo ~ 5. Surprisingly, the static dielectric 
constant remains almost isotropic in the three directions 
of space (c-axis, and ab-plane of cuprates), which means 
that the phonons seem to act in the same manner in these 
three directions. However, it is known that the transport 
properties are incoherent in the c-direction of Cu02 plane. 
A possible way to simulate such a specific situation, is to 
consider a system where the electron bare band mass m* 
is highly anisotropic, while the Coulomb interactions (and 
the electron-phonon interaction) remain isotropic. This 
problem can be studied by considering a crystallized state 
of Frdhlich polarons, where the localizing potential in the 
Wigner approximation is spherically symmetric, but the 



If we consider a polaron mass Mp « lOme and a density electron effective band masses satisfy b = m,j.y/mz ^ 1 



Ts = 100, the resulting Tdeg is of a few Kelvin, while Tc can 
be as high as some fraction of lolq (typically ~ IQQK). 
In conclusion, the phase diagram at strong coupling sep- 
arates into three different regions: polaron crystal (at low 
density, low temperature), "classical" polaron liquid (at 
low density, high temperature), and some liquid state at 
high density where the polarons (or at least some of them) 
are ionized. 

Let us now analyse the situation for a lower than a*, 

where ric < ric (Fig-§, a = 7.2). In that case, by in- 
creasing either the density or the temperature, the crystal 
melting is always driven by criterion (i). Again, the phase 



xy/ 

If we assume & — > 0, which is equivalent to confining 
the carriers in 2D layers, it is known that the polaron 
bound states become anisotropic (m_ this limit, they are 
usually called "surface" polarons Q). Correspondingly, 
the expression (pq) for the variational energy at zero tem- 
perature must be replaced by: 



E = Co + («! + a2 — w) 



a TT 
7^2 



(^ 



2(ai + a2)w'^ 



dt 



° v^Fo" 



o-ait\ 



(1 



=(47) 



o-a2t\ 
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If we compare tliis expression with (Bq) , we see that apart 
from a numerical coefficient, the function to be minimized 
in 2D is equivalent to the one in the 3D case, provided that 
we multiply a by a factor 37r/4 (such a scaling relation 
was first derived by Peeters et al. |32| in the case of an 
isolated polaron). In our case, it is easy to see that, at a 
given density, the variational parameters in 2D are given 
by the simple relations 



2D 



(a) 



,,3-D 



(37ra/4) 



af^(a)=af^(3™/4) 

so that the phase diagram in the anisotropic case can be 
easily deduced from the 3D results (see Fig.0, grey lines 
|]33| ). The main conclusions are summarized in Table 1, 
for two different values of the electron band mass. We 



isotropic (3D) 
anisotropic (2D) 



JJ^ 



(ra* 



jsy 



7.5 
3.2 



6 ■ 10'** 
8 ■ 10^9 



(m* — 2m^) 



5-10' 
6-10^ 



Table 1. Critical density n" {cm ^) and coupling a* (see 
text) Isil. 



first observe that the critical coupling above which po- 
laron dissociation becomes effective is strongly reduced 
from a* (31?) ~ 7.5 (which corresponds to ivlq ~ dtneV 
if we take m* = me), to a* {2D) « 3.2 {ujlo ~ 4:0meV). 
This is reminiscent of the fact that the e-ph coupling has a 
stronger effect in lower dimensions. Moreover, we see that 

the critical density Uc is increased: with the parameters 
of La2Cu04 (m* = 2me, e^ = 30, £00 = 5 pO,^^), one 

obtains ric ~ 6- lQ^'^cm~^, in qualitative agreement with 
experiments. These results suggest that the mechanism of 
polaron dissociation could be relevant in the description 
of the IMT in the high-T^ cuprates. Note however that 
it does not mean that the electron motion would become 
coherent in the c-direction above the dissociation. The re- 
sulting state above the dissociation is out of the goal of 
the present paper. 



5 Conclusion 

In this work, we have attacked the problem of the insulator- 
to-metal transition in doped polar semiconductors. Ne- 
glecting disorder, and replacing the doping ions by a rigid 
compensating jellium, the two main ingredients in the 
problem are: (i) the electron-phonon interaction, which 
leads to polaron formation and (ii) the Coulomb repul- 
sion between the carriers, which causes their crystalliza- 
tion at low density. We have studied the properties of the 
insulating crystallized state in the framework of Wigner 's 
mean-field approximation. The results show that, upon in- 
creasing the Frohlich e-ph coupling a, the Wigner crystal 
of electrons evolves continuously towards a polaron crys- 
tal. On the other hand, a definite distinction between the 



weak and strong e-ph coupling limits comes from the melt- 
ing mechanism occurring for increasing density. By apply- 
ing the Lindemann criterion both to the polaron centcr-of- 
mass R and to the internal coordinate r, we have shown 
that, at T = 0, the crystal melting due to the fluctua- 
tions of the localized particles is only possible for a < a* . 
In the opposite situation a > a* , the IMT is driven by 
polaron dissociation: the particles become more and more 
localized due to the increase in their effective mass, and 
the only way to obtain a metallic state is to ionize the po- 
laron bound states. We have also analysed the melting at 
finite temperature, and we have proposed an approximate 
T vs. n phase diagram for the many-polaron problem. Fi- 
nally, we have shown that if the carriers are confined in 
2D layers, the effect of the e-ph interactions is strongly 
enhanced, so that the mechanism of polaron dissociation 
can be relevant at rather moderate couplings (a ~ 3). 

The overall picture which comes out from this work 
suggests that large polarons should play an important 
role in the superconducting cuprates. Moreover, it will be 
shown in a forthcoming paper |17|| by fully taking into 
account the dipole-dipole interactions between polarons 
that, due to the peculiar dielectric properties of the PWC, 
if some itinerant carriers could coexist with localized po- 
larons, the long-range part of their mutual interactions 
would be overscreened. This would lead to a supercon- 
ducting instability rather than to an insulator-to-metal 
transition, as was already pointed out in reference [ |l5| . 

It is also interesting to compare our results on the melt- 
ing of a PWC to the studies carried out by De Filippis 
et al. for the same model, but in the metallic phase g. 
These authors found that at intermediate e-ph coupling 
(a « 6), the metallic phase can become unstable with re- 
spect to the formation of a Charge Density Wave (CDW), 
as the density is decreased down to n « 10^*~^°c77i~^. For 
the same value of a, we find that in the same range of 
densities, the PWC becomes unstable with respect to the 
formation of a polaron liquid. In this sense, both studies 
on the same problem, either starting from the low density 
limit or from the high density metallic phase, seem to be 
quite consistent. Nevertheless, our analysis suggests that 
the melting of the PWC at zero temperature is first or- 
der, as is the case for the Mott transiton (the free carrier 
density is a discontinuous quantity at the transition). On 
the contrary, the occurence of a CDW instability should 
rather be second order (the amplitude of the component 
Pq of the electronic density with CDW-wavevector q be- 
ing the order parameter). We believe that at stronger cou- 
pling, where the polaron dissociation holds, there should 
exist a metastability of the metallic phase with respect 
to the PWC. This could be related to the fact that the 
long-range Coulomb interactions (i.e. both the electron- 
electron and the electron-phonon interactions) are com- 
pletely unscreened in the crystallized state. 

The behaviour of the long-range Coulomb forces being 
an essential ingredient to the understanding of this "in- 
termediate density" physics is also suggested by a recent 
work by Leggett Q . Arguing, on the basis of experimen- 
tal data on the cuprates, that the Coulomb energy at the 
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metal-superconducting transition is saved in the mid infra- 
red range at small q vectors, Leggett has pointed out that 
the basic pairing mechanism in high-Tc superconductors 
should also involve small in-plane q vectors, and thus could 
be due to long-range interactions. 

Finally, we would like to mention that the present the- 
oretical work can also find applications in different classes 
of physical systems. In particular, it was shown by Jackson 
and Platzman p^ that the dynamics of ripplons coupled 
to a 2D electron set on a film of liquid helium can be 
viewed as a 2D Frohlich polaron problem. 

This work received financial support from the Euro- 
pean Commission (contract no. ERBFMBICT 961230). 



A Details of the Feynman treatment 

Here we evaluate the upper bound to the free energy F 
according to eq. (p3[), that we rewrite here as 



the last expression and put /(r) = ikx5{T^t)~ikx5{T—s), 
eq. (53) takes the form 



F = Fq + A + B (48) 



A.l Evaluation of the variational free energy F 

After the diagonalization has been carried out, the model 
( |22| ) represents two independent harmonic oscillators in 
three dimensions. The corresponding partition function is 
then 

Zq = (2sinh/3ai/2)"3(2sinh/3a2/2)"^ (49) 

This result must be equivalent to the path integral of the 
electronic action (ETI) , including the part corresponding to 
the motion of the auxiliary particle X. In other words. 



Zo = (2sinh/3w/2)"3 / X'(pat/i)e^''["l 



(50) 



If we equate these two results, and make use of the defi- 
nition (p4|), we find: 



^o^^log 



smh(3w/2 



2sinh/3ai/2sinh/3a2/2 



(51) 



I{kx,t,s) ex / 'D{path) exp —— x^dr — — ^ / x^dr 

-^JJ\x{T)-x{a)\^G^{P,T^a)dTd<j 

+ J /(r) • xiT)dT (54) 

The calculation of such a gaussian integral is standard: 
first evaluate the function x{t) for which the exponent 
is maximum, subject to the boundary condition a;(0) = 
x{P) = 0, then (within an unimportant constant) the in- 
tegral (Q) reduces to 



2./ 2e 



'dr 



I{kx,t,s) = exp 
- ^ f fW)-x{a)\^G^{(3,T^a)dTda 



+ J f{T)x{T)dT 

which can be simplified to 



(55) 



/(fc,, t, s) = exp I i y /(r) xiT)dT\ = e"^ [^(*)-*(^)l 
The function x{t) is the solution of the integral equation 



i(r) = 2C / [x{t) - x(cr)]e-'"l^-'^lda -t- -^x{t) - /(r) 
which can be solved by defining the auxiliary function 



X{r) 



w 



G^(/3, T — a)x{a)da 



(57) 



Let us now calculate the term B = (Sq) /I3, which is 
defined as 

^ = -S/ / dtdsG^{p,t-s)(\u{t)-u{s)\^) (52) 
^P Jo Jo 

where G^W, t — s) is the phonon propagator defined in 
section 2.1, and C = Kw/A. The integrand in B can be 
obtained by expanding up to 2nd order in k^ the expres- 
sion 

I(kx,t,s)^ /g.fc.[x(t)-:.(s)]\ (53^ 

{x is one of the three equivalent cartesian coordinates: 
u = (x, ?;,z)). If we write explicitly the path integral in 



Now the differential system reads 

A /^ 2 

x{t) = — W) - X{t)] + ^x(i) - fit) (58) 

W Es 

X{t) ^ v?[X{t) - x{t)] 

In fact, the problem is equivalent to solving the equations 
of motion of the two-body model (|2^) under the external 
driving force /(r). After some lenghty algebra, we find 

/(fc„i,5)=exp|-^.g(i-s)| (59) 



ith 



and 



g{t~s) = —g^{t~s) + ^g2{t~s) (60) 

ai a2 



. . coshl3ai/2 - cosha^it - s ~ 13/2) , 
gi{t-s)^ . , ^ ; 1 = 1,2 



sinh j3ai 



(61) 
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In the isotropic case, by exp anding to second order in k 
both equations (|3|) and (|5|) , with k"^ = kl + ky + kl (the 
three directions are equivalent) and performing the time 
integrations in eq. ( |52| ) we obtain 



B 



iC 



w(a. 



Pct2 Pcti 

a2 coth Qfi coth 

2 2 



The term A = — (S) / [3 is defined as 



A^ 



V8(3Jo 



(|u(t) - u{s)\-^) G^^^{[3,t- s)dtds 



It can be calculated by introducing the Fourier transform 



{\n{t)-n{s)\-') 
(fk 



d^k 



^Jk-[u(t)-u(s)]x 



27r2/c2 ^ 
Iik^,t^s)Iiky,t-s)Iik,,t-s) (65) 



27r2fc2 

The integration over k gives the result 

a 1 [^ dte-* 

A = 



v^i-e-^^o ^fm 



(66) 



If we take the limit /? ^ oo, and include the constant 
part Co = —{^e^a^ /^esrs){me/'m*), the upper bound to 
the ground-state energy takes the form 






a 



e-^dt 



V^Jo ,/|i(l_e-.*) + |j(l-e-^*) 



(67) 



A. 2 Evaluation of the internal and external 
fluctuations 

Here we calculate the average displacements of the center- 
of-mass and of the relative coordinate, which enter in the 
Lindemann criteria (i) and (m). In each space direction, 
since the eigenmodes are independent , we have (^1^2) = 
(ui)(u2) so that making use of eq. (27) and (pO|), we get 
the dimentionlcss quantities 



{SR^ 



{2m*ujLo) ^ 

{2m*u;Lo)~^ 
where 



m 
Mp 



2 r 



^3«+^B« 



at 



al 



a\ 



{5ul) 



{5% 



{5ul) = (2Aiai)-icoth^ 
/3a2 



{5ul) = (2A2a2)"^coth- 



(68) 
(69) 



are the r.m.s. displacements of two harmonic oscillators 
at finite temperature. Observing that the polaronic length 
unit is {2m*ijjLo)^^ — {iTi/m*)eaao, we can write 



(62) 


Rs 

X 

Rs 

X 


V 4a 7 Mp 




./^^f"^)'coth^"^+^^f"M\oth^"^ 
\l ai \w J 2 a2 \w J 2 

/ 2 \ 1/3 1 




U« y al-aj " ^^^> 


') dt ds 
(63) 


>' coth ^"^ + f coth ^"^ 
V ^1 2 v42 2 
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